The Role of (Co)variation in Shaping the Response to Selection in New World Leaf-Nosed Bats.

AbstractUnderstanding and predicting the evolutionary responses of complex morphological traits to selection remains a major challenge in evolutionary biology. Because traits are genetically correlated, selection on a particular trait produces both direct effects on the distribution of that trait and indirect effects on other traits in the population. The correlations between traits can strongly impact evolutionary responses to selection and may thus impose constraints on adaptation. Here, we used museum specimens and comparative quantitative genetic approaches to investigate whether the covariation among cranial traits facilitated or constrained the response to selection during the major dietary transitions in one of the world's most ecologically diverse mammalian families-the phyllostomid bats. We reconstructed the set of net selection gradients that would have acted on each cranial trait during the major transitions to feeding specializations and decomposed the selection responses into their direct and indirect components. We found that for all transitions, most traits capturing craniofacial length evolved toward adaptive directions owing to direct selection. Additionally, we showed instances of dietary transitions in which the complex interaction between the patterns of covariation among traits and the strength and direction of selection either constrained or facilitated evolution. Our work highlights the importance of considering the within-species covariation estimates to quantify evolvability and to disentangle the relative contribution of variational constraints versus selective causes for observed patterns.

Morphological integration during postnatal ontogeny: implications for evolutionary biology.

How covariance patterns of phenotypes change during development is fundamental for a broader understanding of evolution. There is compelling evidence that mammalian cranium covariance patterns change during ontogeny. However, it is unclear to what extent variation in covariance patterns during ontogeny can impact the response to selection. To tackle this question, we explored: (a) the extent to which covariance patterns change during postnatal ontogeny; (b) in which ontogenetic stages covariance patterns differ the most; and (c) the extent to which the phenotypic covariance pattern at different ontogenetic stages can be explained by the same processes determining additive genetic covariance. We sampled the postnatal ontogenetic series for both marsupials and placentals. Within each ontogenetic series, we compared covariance matrices (P-matrices) at different ontogenetic stages. Furthermore, we compared these P-matrices to two target matrices [adult P-matrix and an additive genetic covariance matrix (G-matrix)]. Our results show that for all ontogenetic series, covariance patterns from weaning onward are conserved and probably shaped by the same processes determining the G-matrix. We conclude that irrespective of eventual differences in how selection operates during most of the postnatal ontogeny, the net response to such pressures will probably not be affected by ontogenetic differences in the covariance pattern.

Development and function explain the modular evolution of phalanges in gecko lizards.

Selective regimes favouring the evolution of functional specialization probably affect covariation among phenotypic traits. Phalanges of most tetrapods develop from a conserved module that constrains their relative proportions. In geckos, however, biomechanical specializations associated with adhesive toepads involve morphological variation in the autopodium and might reorganize such modular structures. We tested two hypotheses to explain the modular architecture of hand bones in geckos, one based on developmental interactions and another incorporating functional associations related to locomotion, and compared the empirical support for each hypothetical module between padded and padless lineages. We found strong evidence for developmental modules in most species, which probably reflects embryological constraints during phalangeal formation. Although padded geckos exhibit a functional specialization involving the hyperextension of the distal phalanges that is absent in padless species, the padless species are the ones that show a distal functional module with high integration. Some ancestrally padless geckos apparently deviate from developmental predictions and present a relatively weak developmental module of phalanges and a strongly integrated distal module, which may reflect selective regimes involving incipient frictional adhesion in digit morphology. Modularity of digit elements seems dynamic along the evolutionary history of geckos, being associated with the presence/absence of adhesive toepads.

The pre-eminent role of directional selection in generating extreme morphological change in glyptodonts (Cingulata; Xenarthra).

The prevalence of stasis on macroevolution has been classically taken as evidence of the strong role of stabilizing selection in constraining morphological change. Rates of evolution calculated over longer timescales tend to fall below the expected under genetic drift, suggesting that directional selection signals are erased at longer timescales. Here, we investigated the rates of morphological evolution of the skull in a fossil lineage that underwent extreme morphological modification, the glyptodonts. Contrary to what was expected, we show here that directional selection was the primary process during the evolution of glyptodonts. Furthermore, the reconstruction of selection patterns shows that traits selected to generate a glyptodont morphology are markedly different from those operating on extant armadillos. Changes in both direction and magnitude of selection are probably tied to glyptodonts' invasion of a specialist-herbivore adaptive zone. These results suggest that directional selection might have played a more critical role in the evolution of extreme morphologies than previously imagined.

Detecting patterns of correlational selection with sampling error: A simulation study.

The adoption of a multivariate perspective of selection implies the existence of multivariate adaptive peaks and pervasive correlational selection that promotes co-adaptation between traits. However, to test for the ubiquity of correlational selection in nature, we must first have a sense of how well can we estimate multivariate nonlinear selection (i.e., the γ-matrix) in the face of sampling error. To explore the sampling properties of estimated γ-matrices, we simulated inidividual traits and fitness under a wide range of sample sizes, using different strengths of correlational selection and of stabilizing selection, combined with different number of traits under selection, different amounts of residual variance in fitness, and distinct patterns of selection. We then ran nonlinear regressions with these simulated datasets to simulate γ-matrices after adding random error to individual fitness. To test how well could we detect the imposed pattern of correlational selection at different sample sizes, we measured the similarity between simulated and imposed γ-matrices. We show that detection of the pattern of correlational selection is highly dependent on the total strength of selection on traits and on the amount of residual variance in fitness. Minimum sample size needs to be at least 500 to precisely estimate the pattern of correlational selection. Furthermore, a pattern of selection in which different sets of traits contribute to different functions is the easiest to diagnose, even when using a large number of traits (10 traits), but with sample sizes in the order of 1000 individuals. Consequently, we recommend working with sets of traits from distinct functional complexes and fitness proxies less prone to effects of environmental and demographic stochasticity to test for correlational selection with lower sample sizes.

Measuring the magnitude of morphological integration: The effect of differences in morphometric representations and the inclusion of size.

The magnitude of morphological integration is a major aspect of multivariate evolution, providing a simple measure of the intensity of association between morphological traits. Studies concerned with morphological integration usually translate phenotypes into morphometric representations to quantify how different morphological elements covary. Geometric and classic morphometric representations translate biological form in different ways, raising the question if magnitudes of morphological integration estimates obtained from different morphometric representations are compatible. Here we sought to answer this question using the relative eigenvalue variance of the covariance matrix obtained for both geometric and classical representations of empirical and simulated datasets. We quantified the magnitude of morphological integration for both shape and form and compared results between representations. Furthermore, we compared integration values between shape and form to evaluate the effect of the inclusion or not of size on the quantification of the magnitude of morphological integration. Results show that the choice of morphological representation has significant impact on the integration magnitude estimate, either for shape or form. Despite this, ordination of the integration values within representations is relatively the same, allowing for similar conclusions to be reached using different methods. However, the inclusion of size in the dataset significantly changes the estimates of magnitude of morphological integration, hindering the comparison of this statistic obtained from different spaces. Morphometricians should be aware of these differences and must consider how biological hypothesis translate into predictions about integration in each particular choice of representation.

Genomic Perspective on Multivariate Variation, Pleiotropy, and Evolution.

Multivariate quantitative genetics provides a powerful framework for understanding patterns and processes of phenotypic evolution. Quantitative genetics parameters, like trait heritability or the G-matrix for sets of traits, can be used to predict evolutionary response or to understand the evolutionary history of a population. These population-level approaches have proven to be extremely successful, but the underlying genetics of multivariate variation and evolutionary change typically remain a black box. Establishing a deeper empirical understanding of how individual genetic effects lead to genetic (co)variation is then crucial to our understanding of the evolutionary process. To delve into this black box, we exploit an experimental population of mice composed from lineages derived by artificial selection. We develop an approach to estimate the multivariate effect of loci and characterize these vectors of effects in terms of their magnitude and alignment with the direction of evolutionary divergence. Using these estimates, we reconstruct the traits in the ancestral populations and quantify how much of the divergence is due to genetic effects. Finally, we also use these vectors to decompose patterns of genetic covariation and examine the relationship between these components and the corresponding distribution of pleiotropic effects. We find that additive effects are much larger than dominance effects and are more closely aligned with the direction of selection and divergence, with larger effects being more aligned than smaller effects. Pleiotropic effects are highly variable but are, on average, modular. These results are consistent with pleiotropy being partly shaped by selection while reflecting underlying developmental constraints.

A multiple peak adaptive landscape based on feeding strategies and roosting ecology shaped the evolution of cranial covariance structure and morphological differentiation in phyllostomid bats.

We explored the evolution of morphological integration in the most noteworthy example of adaptive radiation in mammals, the New World leaf-nosed bats, using a massive dataset and by combining phylogenetic comparative methods and quantitative genetic approaches. We demonstrated that the phenotypic covariance structure remained conserved on a broader phylogenetic scale but also showed a substantial divergence between interclade comparisons. Most of the phylogenetic structure in the integration space can be explained by splits at the beginning of the diversification of major clades. Our results provide evidence for a multiple peak adaptive landscape in the evolution of cranial covariance structure and morphological differentiation, based upon diet and roosting ecology. In this scenario, the successful radiation of phyllostomid bats was triggered by the diversification of dietary and roosting strategies, and the invasion of these new adaptive zones lead to changes in phenotypic covariance structure and average morphology. Our results suggest that intense natural selection preceded the invasion of these new adaptive zones and played a fundamental role in shaping cranial covariance structure and morphological differentiation in this hyperdiverse clade of mammals. Finally, our study demonstrates the power of combining comparative methods and quantitative genetic approaches when investigating the evolution of complex morphologies.

Contrasting patterns of RUNX2 repeat variations are associated with palate shape in phyllostomid bats and New World primates.

Establishing the genetic basis that underlies craniofacial variability in natural populations is one of the main topics of evolutionary and developmental studies. One of the genes associated with mammal craniofacial variability is RUNX2, and in the present study we investigated the association between craniofacial length and width and RUNX2 across New World bats (Phyllostomidae) and primates (Catarrhini and Platyrrhini). Our results showed contrasting patterns of association between the glutamate/alanine ratios (Q/A ratio) and palate shape in these highly diverse groups. In phyllostomid bats, we found an association between shorter/broader faces and increase of the Q/A ratio. In New World monkeys (NWM) there was a positive correlation of increasing Q/A ratios to more elongated faces. Our findings reinforced the role of the Q/A ratio as a flexible genetic mechanism that would rapidly change the time of skull ossification throughout development. However, we propose a scenario in which the influence of this genetic adjustment system is indirect. The Q/A ratio would not lead to a specific phenotype, but throughout the history of a lineage, would act along with evolutionary constraints, as well as other genes, as a facilitator for adaptive morphological changes.

Evolution of morphological integration in the skull of Carnivora (Mammalia): Changes in Canidae lead to increased evolutionary potential of facial traits.

Morphological integration refers to the fact that different phenotypic traits of organisms are not fully independent from each other, and tend to covary to different degrees. The covariation among traits is thought to reflect properties of the species' genetic architecture and thus can have an impact on evolutionary responses. Furthermore, if morphological integration changes along the history of a group, inferences of past selection regimes might be problematic. Here, we evaluated the stability and evolution of the morphological integration of skull traits in Carnivora by using evolutionary simulations and phylogenetic comparative methods. Our results show that carnivoran species are able to respond to natural selection in a very similar way. Our comparative analyses show that the phylogenetic signal for pattern of integration is lower than that observed for morphology (trait averages), and that integration was stable throughout the evolution of the group. That notwithstanding, Canidae differed from other families by having higher integration, evolvability, flexibility, and allometric coefficients on the facial region. These changes might have allowed canids to rapidly adapt to different food sources, helping to explain not only the phenotypic diversification of the family, but also why humans were able to generate such a great diversity of dog breeds through artificial selection.

Intense natural selection preceded the invasion of new adaptive zones during the radiation of New World leaf-nosed bats.

The family Phyllostomidae, which evolved in the New World during the last 30 million years, represents one of the largest and most morphologically diverse mammal families. Due to its uniquely diverse functional morphology, the phyllostomid skull is presumed to have evolved under strong directional selection; however, quantitative estimation of the strength of selection in this extraordinary lineage has not been reported. Here, we used comparative quantitative genetics approaches to elucidate the processes that drove cranial evolution in phyllostomids. We also quantified the strength of selection and explored its association with dietary transitions and specialization along the phyllostomid phylogeny. Our results suggest that natural selection was the evolutionary process responsible for cranial diversification in phyllostomid bats. Remarkably, the strongest selection in the phyllostomid phylogeny was associated with dietary specialization and the origination of novel feeding habits, suggesting that the adaptive diversification of phyllostomid bats was triggered by ecological opportunities. These findings are consistent with Simpson's quantum evolutionary model of transitions between adaptive zones. The multivariate analyses used in this study provides a powerful tool for understanding the role of evolutionary processes in shaping phenotypic diversity in any group on both micro- and macroevolutionary scales.

The evolution of phenotypic integration: How directional selection reshapes covariation in mice.

Variation is the basis for evolution, and understanding how variation can evolve is a central question in biology. In complex phenotypes, covariation plays an even more important role, as genetic associations between traits can bias and alter evolutionary change. Covariation can be shaped by complex interactions between loci, and this genetic architecture can also change during evolution. In this article, we analyzed mouse lines experimentally selected for changes in size to address the question of how multivariate covariation changes under directional selection, as well as to identify the consequences of these changes to evolution. Selected lines showed a clear restructuring of covariation in their cranium and, instead of depleting their size variation, these lines increased their magnitude of integration and the proportion of variation associated with the direction of selection. This result is compatible with recent theoretical works on the evolution of covariation that take the complexities of genetic architecture into account. This result also contradicts the traditional view of the effects of selection on available covariation and suggests a much more complex view of how populations respond to selection.

Timing and patterns of diversification in the Neotropical bat genus Pteronotus (Mormoopidae).

We investigate the biogeographic processes related to the origin and current patterns of distribution of the extant species of the genus Pteronotus. This clade of insectivorous bats is widely distributed in the Neotropical Region and has recently gone through a taxonomic update which increased more than twice its diversity. Using six molecular markers of 15 Pteronotus lineages ranging from Mexico to Central Brazil, we reconstruct a time-calibrated tree and infer molecular evolutionary rates for this bat genus. In addition, estimates of range evolution across phylogeny were obtained through statistical model testing among six different biogeographic models. The origin of the genus Pteronotus occurred approximately 16million years ago (Ma), with initial cladogenesis events being evenly distributed across the phylogeny. Divergence between most closely related species is recent, falling in the Pleistocene period less than 2.6Ma. Mainland lineages present congruent patterns of north versus south continent splitting while insular clades differ in their time of arrival in the Caribbean Islands. Temporal and geographic range estimates for early nodes of Pteronotus phylogeny suggest a central role of Neogene tectonic reorganizations of Central America in the group diversification process. Also, South American colonization by Pteronotus occurred early in the genus history. Founder-event speciation was an important mode of lineage splitting in Pteronotus, with two independent dispersal jumps having occurred to the Greater Antilles. Finally, Pleistocenic sea-level variation and climatic oscillations are possibly associated with divergence between sister-species and recent ages of MRCA for Pteronotus species.

Evolutionary processes and its environmental correlates in the cranial morphology of western chipmunks (Tamias).

The importance of the environment in shaping phenotypic evolution lies at the core of evolutionary biology. Chipmunks of the genus Tamias (subgenus Neotamias) are part of a very recent radiation, occupying a wide range of environments with marked niche partitioning among species. One open question is if and how those differences in environments affected phenotypic evolution in this lineage. Herein we examine the relative importance of genetic drift versus natural selection in the origin of cranial diversity exhibited by clade members. We also explore the degree to which variation in potential selective agents (environmental variables) are correlated with the patterns of morphological variation presented. We found that genetic drift cannot explain morphological diversification in the group, thus supporting the potential role of natural selection as the predominant evolutionary force during Neotamias cranial diversification, although the strength of selection varied greatly among species. This morphological diversification, in turn, was correlated with environmental conditions, suggesting a possible causal relationship. These results underscore that extant Neotamias represent a radiation in which aspects of the environment might have acted as the selective force driving species' divergence.

Evolution of a complex phenotype with biphasic ontogeny: Contribution of development versus function and climatic variation to skull modularity in toads.

The theory of morphological integration and modularity predicts that if functional correlations among traits are relevant to mean population fitness, the genetic basis of development will be molded by stabilizing selection to match functional patterns. Yet, how much functional interactions actually shape the fitness landscape is still an open question. We used the anuran skull as a model of a complex phenotype for which we can separate developmental and functional modularity. We hypothesized that functional modularity associated to functional demands of the adult skull would overcome developmental modularity associated to bone origin at the larval phase because metamorphosis would erase the developmental signal. We tested this hypothesis in toad species of the Rhinella granulosa complex using species phenotypic correlation pattern (P-matrices). Given that the toad species are distributed in very distinct habitats and the skull has important functions related to climatic conditions, we also hypothesized that differences in skull trait covariance pattern are associated to differences in climatic variables among species. Functional and hormonal-regulated modules are more conspicuous than developmental modules only when size variation is retained on species P-matrices. Without size variation, there is a clear modularity signal of developmental units, but most species have the functional model as the best supported by empirical data without allometric size variation. Closely related toad species have more similar climatic niches and P-matrices than distantly related species, suggesting phylogenetic niche conservatism. We infer that the modularity signal due to embryonic origin of bones, which happens early in ontogeny, is blurred by the process of growth that occurs later in ontogeny. We suggest that the species differing in the preferred modularity model have different demands on the orbital functional unit and that species contrasting in climate are subjected to divergent patterns of natural selection associated to neurocranial allometry and T3 hormone regulation.

High evolutionary constraints limited adaptive responses to past climate changes in toad skulls.

Interactions among traits that build a complex structure may be represented as genetic covariation and correlation. Genetic correlations may act as constraints, deflecting the evolutionary response from the direction of natural selection. We investigated the relative importance of drift, selection, and constraints in driving skull divergence in a group of related toad species. The distributional range of these species encompasses very distinct habitats with important climatic differences and the species are primarily distinguished by differences in their skulls. Some parts of the toad skull, such as the snout, may have functional relevance in reproductive ecology, detecting water cues. Thus, we hypothesized that the species skull divergence was driven by natural selection associated with climatic variation. However, given that all species present high correlations among skull traits, our second prediction was of high constraints deflecting the response to selection. We first extracted the main morphological direction that is expected to be subjected to selection by using within- and between-species covariance matrices. We then used evolutionary regressions to investigate whether divergence along this direction is explained by climatic variation between species. We also used quantitative genetics models to test for a role of random drift versus natural selection in skull divergence and to reconstruct selection gradients along species phylogeny. Climatic variables explained high proportions of between-species variation in the most selected axis. However, most evolutionary responses were not in the direction of selection, but aligned with the direction of allometric size, the dimension of highest phenotypic variance in the ancestral population. We conclude that toad species have responded to selection related to climate in their skulls, yet high evolutionary constraints dominated species divergence and may limit species responses to future climate change.

Evolution of the Genotype-to-Phenotype Map and the Cost of Pleiotropy in Mammals.

Evolutionary studies have long emphasized that the genetic architecture of traits holds important microevolutionary consequences. Yet, studies comparing the genetic architecture of traits across species are rare, and discussions of the evolution of genetic systems are made on theoretical arguments rather than on empirical evidence. Here, we compared the genetic architecture of cranial traits in two different mammalian model organisms: the gray short-tailed opossum, Monodelphis domestica, and the laboratory mouse, Mus musculus We show that both organisms share a highly polygenic genetic architecture for craniofacial traits, with many loci of small effect. However, these two model species differ significantly in the overall degree of pleiotropy, N, of the genotype-to-phenotype map, with opossums presenting a higher average N They also diverge in their degree of genetic modularity, with opossums presenting less modular patterns of genetic association among traits. We argue that such differences highlight the context dependency of gene effects, with developmental systems shaping the variational properties of genetic systems. Finally, we also demonstrate based on the opossum data that current measurements for the relationship between the mutational effect size and N need to be re-evaluated in relation to the importance of the cost of pleiotropy for mammals.

Integrating multiple evidences in taxonomy: species diversity and phylogeny of mustached bats (Mormoopidae: Pteronotus).

A phylogenetic systematic perspective is instrumental in recovering new species and their evolutionary relationships. The advent of new technologies for molecular and morphological data acquisition and analysis, allied to the integration of knowledge from different areas, such as ecology and population genetics, allows for the emergence of more rigorous, accurate and complete scientific hypothesis on species diversity. Mustached bats (genus Pteronotus) are a good model for the application of this integrative approach. They are a widely distributed and a morphologically homogeneous group, but comprising species with remarkable differences in their echolocation strategy and feeding behavior. The latest systematic review suggested six species with 17 subspecies in Pteronotus. Subsequent studies using discrete morphological characters supported the same arrangement. However, recent papers reported high levels of genetic divergence among conspecific taxa followed by bioacoustic and geographic agreement, suggesting an underestimated diversity in the genus. To date, no study merging genetic evidences and morphometric variation along the entire geographic range of this group has been attempted. Based on a comprehensive sampling including representatives of all current taxonomic units, we attempt to delimit species in Pteronotus through the application of multiple methodologies and hierarchically distinct datasets. The molecular approach includes six molecular markers from three genetic transmission systems; morphological investigations used 41 euclidean distances estimated through three-dimensional landmarks collected from 1628 skulls. The phylogenetic analysis reveals a greater diversity than previously reported, with a high correspondence among the genetic lineages and the currently recognized subspecies in the genus. Discriminant analysis of variables describing size and shape of cranial bones support the rising of the genetic groups to the specific status. Based on multiples evidences, we present an updated taxonomic arrangement composed by 16 extant species and a new and more robust phylogenetic hypothesis for the species included in the genus Pteronotus. Studies developed under such integrative taxonomic approach are timely for a deeper and wider comprehension of Neotropical diversity, representing the first step for answering broader questions on evolutionary and ecological aspects of Neotropical life history.

Morphological and dietary responses of chipmunks to a century of climate change.

Predicting how individual taxa will respond to climatic change is challenging, in part because the impacts of environmental conditions can vary markedly, even among closely related species. Studies of chipmunks (Tamias spp.) in Yosemite National Park provide an important opportunity to explore the reasons for this variation in response. While the alpine chipmunk (T. alpinus) has undergone a significant elevational range contraction over the past century, the congeneric and partially sympatric lodgepole chipmunk (T. speciosus) has not experienced an elevational range shift during this period. As a first step toward identifying the factors underlying this difference in response, we examined evidence for dietary changes and changes in cranial morphology in these species over the past century. Stable isotope analyses of fur samples from modern and historical museum specimens of these species collected at the same localities indicated that signatures of dietary change were more pronounced in T. alpinus, although diet breadth did not differ consistently between the study species. Morphometric analyses of crania from these specimens revealed significant changes in cranial shape for T. alpinus, with less pronounced changes in shape for T. speciosus; evidence of selection on skull morphology was detected for T. alpinus, but not for T. speciosus. These results are consistent with growing evidence that T. alpinus is generally more responsive to environmental change than T. speciosus, but emphasize the complex and often geographically variable nature of such responses. Accordingly, future studies that make use of the taxonomically and spatially integrative approach employed here may prove particularly informative regarding relationships between environmental conditions, range changes, and patterns of phenotypic variation.